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GROUND STATE ON THE DUMBBELL GRAPH 


JEREMY L. MARZUOLA AND DMITRY E. PELINOVSKY 


Abstract. We consider standing waves in the focusing nonlinear Schrodinger (NLS) equation on 
a dumbbell graph (two rings attached to a central line segment subject to the Kirchhoff boundary 
conditions at the junctions). In the limit of small norm, the ground state (the orbitally stable 
standing wave of the smallest energy at a fixed norm) is represented by a constant solution. 
However, when the norm is increased, this constant solution undertakes two bifurcations, where 
the first is the pitchfork (symmetry breaking) bifurcation and the second one is the symmetry 
preserving bifurcation. As a result of the first symmetry breaking bifurcation, the standing wave 
becomes more localized in one of the two rings. As a result of the second symmetry preserving 
bifurcation, the standing wave becomes localized in the central line segment. In the limit of large 
norm solutions, both standing waves are represented by a truncated solitary wave localized in 
either the ring or the central line segment. Although the asymmetric wave supported in the ring is 
a ground state near the symmetry breaking bifurcation of the constant solution, it is the symmetric 
wave supported in the central line segment which becomes the ground state in the limit of large 
norm. The analytical results are confirmed by numerical approximations of the ground state on 
the dumbbell graph. 


1. Introduction 

Nonlinear Schrodinger (NLS) equations on quantum graphs have been recently studied in many 
physical and mathematical aspects [15]. In the physical literature, mostly in the context of Bose- 
Einstein condensation, various types of graphs have been modeled to show formation and trapping of 
standing waves [61II1I2511261130]. In the mathematical literature, existence, variational properties, 
stability, and scattering have been studied on star graphs, including the U-shaped graphs Enam. 

More complicated graphs may lead to resonances and nontrivial bifurcations of standing waves. 
For example, standing waves were studied on the tadpole graph (a ring attached to a semi-inhnite 
line) [71[T9]. Besides the standing waves supported in the ring that bifurcate from eigenvalues of the 
linear operators closed in the ring, the tadpole graph also admits the standing waves localized in 
the ring with the tails extended in the semi-inhnite line. These standing waves bifurcate from the 
end-point resonance of the linear operators dehned on the tadpole graph and include the positive 
solution, which is proved to be orbitally stable in the evolution of the cubic NLS equation near the 
bifurcation point m- The positive solution bifurcating from the end-point resonance bears the 
lowest energy at the hxed norm, called the ground state. Other positive states on the tadpole 
graph also exist in parameter space far away from the end-point resonance but they do not bear 
smallest energy and they do not branch off the ground state [7j- 

The present contribution is devoted to analysis of standing waves (and the ground state) on the 
dumbbell graph shown in Figure [TJ The dumbbell graph is represented by two rings (with equal 
length normalized to 27r) connected by the central line segment (with length 2L). At the junctions 
between the rings and the central line segment, we apply the Kirchhoff boundary conditions to 
dehne the coupling. These boundary conditions ensure continuity of functions and conservation 
of the current flow through the network junctions; they also allow for self-adjoint extension of the 
Laplacian operator dehned on the dumbbell graph. 
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Figure 1 . Schematic representation of the dumbbell graph. 


Let the central line segment 
/_ := [—L — 2tt, —L] and /+ := 


be placed on Iq := [—L,L], whereas the end rings are placed on 
[L, L + 27r]. The Laplacian operator is defined piecewise by 


AT 


n" (x), X G 
Uo(x), X G /o, 
ix'l(x), X G /+, _ 



n_(x) 

X G 

acting on T = 

'Uo(x) 

X G Iq, 


. u+{x) 

X G /+, 


subject to the Kirchhoff boundary conditions at the two junctions: 


( 1 . 1 ) 

and 

( 1 . 2 ) 


U-{—L — 2tt) = U-{—L) = uo{—L), 
u'_{—L) — u'_{—L — 27r) = Uq(—L), 


f u+(L + 27r) = u+(L) = uo(L), 
I tt+(L) - rt+(L + 27r) = Uq(L). 


The Laplacian operator A is equipped with the domain T’(A) given by a subspace of U Iq U 

/+) closed with the boundary conditions (11.11) and (11.21) . By Theorem 1.4.4 in [5], the Kirchhoff 
boundary conditions are symmetric and the operator A is self-adjoint on its domain 'D{A). 

The cubic NLS equation on the dumbbell graph is given by 

d 

(1.3) z—T = AT + 2|T|2 t , TgP(A), 

where the nonlinear term |TpT is also defined piecewise on /_ U Iq U /+. The energy of the cubic 
NLS equation (11.31) is given by 


(1.4) F(T)= [ (|a,T|2-|T|4)dx, 

2/_u/oU/+ 

and it is conserved in the time evolution of the NLS equation (11.31) . The energy is dehned in the 
energy space T(A) given by 


T(A) := 


|t G H\I.UIoUI+) : 


U-{—L — 27r) = u-{—L) = uq{—L) 
tt+(L -|- 27r) = rt+(L) = uq{L) 


Local and global wellposedness of the cubic NLS equation (|1.3I) both in energy space T(A) and 
domain space 11(A) can be proved using standard techniques, see [3]. 
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Standing waves of the focusing NLS equation (11.31) are given by the solutions of the form 'h(t, x) = 
where A and G ’P(A) are considered to be real. This pair satisfies the stationary NLS 

equation 

(1.5) - - 2|^>p^> = A€M,^>GP(A). 

The stationary NLS equation (jl.5l) is the Euler-Lagrange equation of the energy functional := 
E — AQ, where the charge 

(1.6) Q(^') = [ \^\‘^dx 

J1-UI0UI+ 

is another conserved quantity in the time evolution of the NLS equation (|1.3I) . 

We shall now define the ground state of the NLS equation on the dumbbell graph as the standing 
wave of smallest energy E at a fixed value of Q, that is, a solution of the constrained minimization 
problem 

(1.7) Eo = inf{E(^): ^ G f(A), Q(^-) = Qo}- 

By Theorem 1.4.11 in [5], although the energy space f(A) is only defined by the continuity 
boundary conditions, the Kirchhoff boundary conditions for the derivatives are natural boundary 
conditions for critical points of the energy functional E(4') in the space f(A). In other words, 
using test functions and the weak formulation of the Euler-Lagrange equations in the energy space 
£^(A), the derivative boundary conditions are also obtained in addition to the continuity boundary 
conditions. By bootstrapping arguments, we conclude that any critical point of the energy func¬ 
tional H\ in f(A) is also a solution of the stationary NLS equation (II.5p in 'D{A). On the other 
hand, solutions of the stationary NLS equation (II.5p in T>{A) are immediately the critical points 
of the energy functional H\. Therefore, the set of standing wave solutions of the stationary NLS 
equation (jl.5|) is equivalent to the set of critical points of the energy functional H\. 

Non-existence of ground states on graphs was proved by Adami et al. [1] for some non-compact 
graphs. For example, the graph consisting of one ring connected to two semi-infinite lines does not 
have a ground state. On the other hand, the tadpole graph with one ring and one semi-infinite line 
escapes the non-existence condition of [4] and has a ground state, in agreement with the results of 
|19j . Because /_ U Iq U /+ is compact, existence of the global constrained minimizer in (II.7h follows 
from the standard results in calculus of variations. As a minimizer of the energy functional 
the ground state is orbitally stable in the time evolution of the NLS equation (11.311 . see for instance 
m- The main question we would like to answer is how the ground state looks like on the dumbbell 
graph depending on the parameter Qq for the charge Q. Until now, no rigorous analysis of the 
NLS equation (11.311 on a compact graph has been developed. On the other hand, ground states on 
compact intervals subject to Dirichlet or periodic boundary conditions have been considered in the 
literature UM- 

The dumbbell graph resembles the geometric configuration that arises typically in the double¬ 
well potential modeled by the Gross-Pitaevskii equation [ini dH [Ml [Ml [23]. From this analogy, 
one can anticipate that the ground state is a symmetric state distributed everywhere in the graph 
in the limit of small values of Qq but it may become an asymmetric standing wave residing in just 
one ring as a result of a pitchfork bifurcation for larger values of Qq. We show in this paper that 
this intuitive picture is only partially correct. 

We show that the ground state is indeed represented by a constant (symmetric) solution for 
small values of Qq. For larger values of Qo, the constant solution undertakes two instability bi¬ 
furcations. At the first bifurcation associated with the anti-symmetric perturbation, a family of 
positive asymmetric standing waves is generated. The asymmetric wave has the lowest energy at 
the fixed Qo near the symmetry breaking bifurcation. At the second bifurcation associated with 
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the symmetric perturbation of the constant solution, another family of positive symmetric standing 
waves is generated. The symmetric wave does not have the lowest energy at the fixed Qo near the 
bifurcation but has this property in the limit of large Qq. It is rather surprising that both the 
precedence of the symmetry-breaking bifurcation of the constant solution and the appearance of 
the symmetric wave as a ground state in the limit of large Qo do not depend on the value of the 
length parameter L relative to vr. 

Our main result is formulated as the following two theorems. We also include numerical approx¬ 
imations of the standing waves of the stationary NLS equation (11.511 in order to illustrate the main 
result. The numerical work relies on the Petviashvili’s and Newton’s iterative methods which are 
commonly used for approximation of standing waves of the NLS equations [2ll[29]. 

Theorem 1.1. There exist Qq and Qq* ordered as 0 < Qq < Qq* < oo sueh that the ground state 
of the constrained minimization problem ra for Qo G (0, Qq) is given (up to an arbitrary rotation 
of phase) by the constant solution of the stationary NLS equation 

(1.8) ^{x) = p, A = -2p^, Qo = 2{L + 2'ii)p^. 

The constant solution undertakes the symmetry breaking bifurcation at Qq and the symmetry pre¬ 
serving bifureation at Q”, which result in the appearance of new positive non-eonstant solutions. 
The asymmetrie standing wave is a ground state of (Q for Q ^ Qo ^rit the symmetric standing 
wave is not a ground state of for Q > Qo*- 

Theorem 1.2. In the limit of large negative A, there exist two standing wave solutions of the 
stationary NLS equation 11.5|) . One solution is a positive asymmetrie wave localized in the ring: 

(1.9) <l>(x) = |A|^/^sech(|A|^/^(a: - L - vr))-k $(x), Qo = 2|A|^/^-k Qo, 

and the other solution is a positive symmetric wave loealized in the central line segment: 

(1.10) <h(x) = |A|^/^sech(|A|^/^a:)-k i(x), Qo = 2|A|^/^-k Qo, 

where ||^||_ff2(/_u7oU/+) 0 ond |Qo| —5-0 as A —>• —oo in both cases. The positive symmetric wave 

satisfuina di.iOl) is a around state of the constrained minimization vroblem ra for Qo sufficiently 
large. 

Remark 1.3. It follows from Lemmas 13.21 15^ and Remark 13.51 that the constant standing wave 
(|1.8p undertakes a sequence of bifurcations, where the first two bifurcations at Qq and Qq* lead 
to the positive asymmetric and symmetric standing waves respectively. We also show numerically 
that these same positive waves are connected to the truncated solitary waves (|1.9p and (ll.lOp as 
A —)• —oo. See Figures [6l [71 [HI and[9l 

Remark 1.4. We show in Lemma 14.91 that the energy difference between the truncated solitary 
wave localized in the central segment and the one localized in the ring is exponentially small 
as A —?■ —oo. Nevertheless, it is only the former wave which appears to be the ground state 
of the constrained minimization problem (II.7p as Qo —)• oo. In the numerical iterations of the 
Petviashvili’s and Newton’s methods, both standing waves arise naturally when the initial data is 
concentrated either in a loop or in a central link. Figures [TOl and [11] illustrate that although the 
positive asymmetric wave is not a global constrained minimizer of energy, it is a local constrained 
minimizer, so that it is also orbitally stable in the time evolution of the cubic NLS equation ()1.3p . 

The paper is organized as follows. Section 2 reports a complete characterization of the linear 
spectrum of the Laplacian operator on the dumbbell graph. Section 3 is devoted to the analytical 
characterization of the constant standing wave (II.8p and the first two instability bifurcations when 
parameter Qq is increased. Section 4 describes the analytical characterization of the two standing 
waves localized in the central segment and at one of the two rings in the limit of large values of 
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Qo- The proofs of Theorems II.11 and II.21 are furnished by the individual results of Sections 3 and 4. 
Section 5 reports numerical approximations of the standing waves of the stationary NTS equation 

m- 


2. Linear spectrum of the Laplacian on the dumbbell graph 


The linear spectrum of the Laplacian on the dumbbell graph is defined by solutions of the spectral 
problem 

(2.1) -AU = XU, AgM, UeV{A). 

Because /_ U Iq U /+ is compact, the spectrum of —A is purely discrete. Let us denote it by cj(—A). 
Because A is self-adjoint with the domain T’(A) in L^(I_ U Iq U /+), the spectrum a{—A) consists 
of real positive eigenvalues of equal algebraic and geometric multiplicities. The distribution of 
eigenvalues is given by the following result. 

Proposition 2.1. cr(—A) consists of a simple zero eigenvalue with the constant eigenfunction and 
the union of the following three countable sequences of eigenvalues: 

• A sequence of double eigenvalues {n^jngM- The corresponding eigenfunctions are compactly 
supported on either /_ or /+ and are odd with respect to the middle point in I±. 

• A sequence of simple eigenvalues where cOn is given by a positive root of the 

transcendental equation 

(2.2) := 2tan(a;7r) -|- tan(a;L) = 0. 

The corresponding eigenfunctions are distributed everywhere in I_ U Iq U /+ and are even 
with respect to the middle point in Iq. 

• A sequenee of simple eigenvalues where Qn is given by a positive root of the 

transcendental equation 

(2.3) := 2tan(a;7r) — cot{ujL) = 0. 

The eorresponding eigenfunetions are distributed everywhere in /_ U Iq U and are odd 
with respect to the middle point in Iq. 


Proof. Let us decompose U in the components {u-,uq,u+} defined on {/_, Iq, 1+} respectively. We 
first observe the following reduction of the spectral problem (|2.ip : if uq is identically zero, then ri_|_ 
and U- are uncoupled, and each satisfies the over-determined boundary-value problem 


, . J -u±(x) = Au±(x), X G /±, 

^ ^ ^ I «±(^) e 

where II‘^^^ q{I±) denotes the subspace of Il'^^j.{I±) subject to the additional Dirichlet boundary 
conditions at the end points of I±. The over-determined problem (12.4|] can be solved in the space 
of functions which are odd with respect to the middle point in I±. In this way, a complete set of 
solutions of the boundary-value problem (12.41) is given by the set of double eigenvalues 
with two linearly independent eigenfunctions 


(2.5) tt+(x) = sin[n(x — L — tt)], tt_(x) = 0 

and 


(2.6) rt+(x) = 0, u_(x) = sin[n(x-|-L-|-vr)]. 

We next consider other solutions of the spectral problem (12.ip . for which uq is not identically zero. 
By the parity symmetry on /_U/oU/+, the eigenfunctions are either even or odd with respect to the 
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middle point in Iq. Since (t(—A) consists of real positive eigenvalues of equal algebraic and geometric 
multiplicities, we parameterize X = uP. For even functions, we normalize the eigenfunction U by 

(2.7) uo{x) = cos{ujx), X £ Iq, U-{—x) = u+{x), X G/+. 

The most general solution of the differential equation (I2.ip on I_|_ is given by 

(2.8) u+{x) = A cos[cu(x — L — tt)] + B sin[a;(x — L — vr)], 

where the coefficients A and B, as well as the spectral parameter w is found from the Kirchhoff 
boundary conditions in (jl.2ll . However, since u+(T + 27r) = u+(L), we have B ^ 0 if and only if 
sin(7rL<;) = 0. This condition is satisfied for cu = n G N, when the uncoupled eigenfunction (|2.5I1 arise 
for the H-term of the decomposition ()2.8p . Therefore, without loss of generality, we can consider 
other eigenfunctions by setting B = 0. Then, the Kirchhoff boundary conditions (11.21) yield the 
constraints 

J Hcos(a;7r) = cos(t(jL), 

P 2Ha; sin(t<;7r) = —a; sin(a;L). 

Eliminating A := we obtain the dispersion relation that admits a countable set 

of positive roots {wnjneN in addition to the zero root cu = 0 that corresponds to the constant 
eigenfunction. 

For odd functions, we normalize the eigenfunction U by 

(2.9) uo(x) = sin(a;x), x G Iq, m_(—x) = —u+(x), x G/+. 

Representing the most general solution of the differential equation (|2.1I) on I_|_ by (|2.8p , we have the 
same reasoning to set B = 0. Then, the Kirchhoff boundary conditions (II.2p yield the constraints 

J Hcos(a;7r) = sin(a;L), 

P 2Auj sin(a;7r) = ui cos{ujL). 

Eliminating A := we obtain the dispersion relation (j2.3l) that admits a countable set of 

positive roots The root cj = 0 is trivial since it corresponds to the zero solution for U. 

Therefore, A = 0 is a simple eigenvalue of the spectral problem (12.ip with constant eigenfunction 
U. All assertions of the proposition are proved. □ 


Remark 2.2. When L is a rational multiplier of one double eigenvalue in the sequence {n^j^eN 
is actually a triple eigenvalue. Indeed, if L = then in addition to the eigenfunctions (12.5p and 
dST]), we obtain the third eigenfunction for X = n?. If m is even, the third eigenfunction is given 

by 



( uo{x) = cos(nx). 

X G Jo, 

(2.10) 

< u+(x) = (—l)"'”''^ cos[n(x — L — tt)]. 

X G /+, 


[ U-{x) = u+(-x). 

X G /- 

whereas if m is odd, the eigenfunction is given by 



no(x) = sin(nx). 

X G Iq, 

(2.11) 

u+(x) = (—I)"'"' ^ cos[n(x — L — vr)]. 

X G /+, 


^ u-{x) = -u+(-x). 

X G /_ 

We say that a resonance 

occurs if L is a rational multiplier of 



Figure [2] show graphical solutions of the dispersion relations (12.2p and (12.3p for T = f with the 
roots {wnIngN and clearly marked. The graphical solution persists for any L < tt as seen 
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Figure 2. Graphical solutions of the dispersion relations (|2.2I1 (top) and (12.311 (bot¬ 
tom) for L = ^. 


from the vertical asymptotics of the functions tan(a;7r), tan(a;L), and cot(a;L). As a result, for 
every L < tt the first positive roots of the dispersion relations (12.2h and (j2.3h satisfy 


( 2 . 12 ) 


0 < rii < - < cji < min 
2 



< 112 < . . ■ 


With similar analysis, it follows that for every L > vr, the first roots satisfy 


(2.13) 


TT 


0 < < — <ui< min 

ZiIj 



< ^2 < • • • 


In either case, the smallest positive eigenvalue flf G cj(—A) corresponds to the odd eigenfunction 
with respect to the middle point in Iq, whereas the second positive eigenvalue G cr(—A) corre¬ 
sponds to the even eigenfunction. As is shown in Section O the order of eigenvalues in (|2.12ll or 
()2.13p is important for analysis of the first two bifurcations of the constant standing wave of the 
stationary NLS equation (II.5p . 

Remark 2.3. Since < ^, < wi < fl 2 , and oji < 1, the first two positive eigenvalues in cr(—A) 

are simple, according to Proposition 12.11 
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For further references, we give the explicit expression for the eigenfunction U corresponding to 
the smallest positive eigenvalue in cr(—A): 


(2.14) 


U{x) = < 


- rs(Sit) cos(f2i(x + L + tt)), xG/_, 

sin(nix), X G Jo, 

— T. — -ttW X G/+. 


^2i^cos(Oi(x-L-7r)), 

V cos(iZi7r) ^ 

The value fli is found from the first positive root of the transcendental equation 
(2.15) 2 sin(cu7r) sin(a;L) = cos(u;7r) cos(wL). 

The following proposition summarizes on properties of the root fli. 


Proposition 2.4. For every L > 0, the first positive root fli of the dispersion relation \2.15\} 
satisfies fli < ^ and fli < ■^. Moreover, fli —^ as L —>■ 0 and 2LVti —>■ vr as L ^ oo. 

Proof. The bounds fli < ^ and fli < ^ follow from the orderings (12.121) and (I2.13p . 

From the bound fli < | and the algebraic equation (j2.15p . we realize that sin(Lfli) —)• 0 and 
cos(7rni) —)■ 0 as L —)• 0. Therefore, fli —^ as L —>• 0. 

On the other hand, from the bound Oi < ^ and the algebraic equation (12.151) . we realize that 
sin(7rni) —)• 0 and cos(LOi) —>■ 0 as L —)• oo. Therefore, 2LOi —>■ tt as L —)• oo. □ 


3. Proof of Theorem 11.11 

As we discussed in the introduction, critical points of the energy fnnctional H\ := E—AQ in £i(A) 
are equivalent to the strong solutions of the stationary NFS equation (|1.5p . Among all solutions 
of the stationary NFS equation (11.51) . we are interested in the solutions that yield the minimum of 
energy denoted by Eq at the fixed charge denoted by Qq. These solutions correspond to the ground 
states of the constrained minimization problem (|1.7p . There exists a mapping between parameters 
A and Qo for the ground state solutions. 

The first elementary result describes the local bifurcation of the constant standing wave (jl.Sp as 
a ground state of the minimization problem (11.71) for small charge Qo- 

Femma 3.1. There exists Qq such that the constant standing wave is a ground state of the 

constrained minimization problem ( 13 ) /or goG (o,gs)- 

Proof. The spectrum of —A equipped with the domain P(A) C L'^{I- U Iq U /+) is given by 
eigenvalues described in Proposition 12.11 The simple zero eigenvalue is the lowest eigenvalue in 
(t(—A) corresponding to the constant eigenfunction. 

By the local bifurcation theory [13 E], the ground state of the constrained minimization problem 
(II.7p for small values of Qq is the standing wave bifurcating from the constant eigenfunction of —A 
in P(A). Since the constant solution $(x) = p exists for the stationary NFS equation (|1.5p with 
A = —2p‘^ for every A < 0, the bifurcating ground state is the constant standing wave given by 
(USD. The relation between go and A is computed explicitly from the definition (|1.5I) : 

go = 2(F + 2'K^pfi = (L + 27r)|A|. 

This concludes the proof of the lemma. □ 

We shall now consider bifurcations of new standing waves of the stationary NFS equation (11.51) 
from the constant solution (|1.8p . For both L < tt and L > vr, the lowest nonzero eigenvalue in 
(t(—A) in Proposition 12.II is the positive eigenvalue 11^, which corresponds to the odd eigenfunction 
in /_ U Jo U /+, see orderings (I2.12p and (I2.13p . This smallest nonzero eigenvalue induces the 
symmetry-breaking bifurcation of the ground state of the constrained minimization problem (11.71) 
at gQ. This bifurcation only marks the first bifurcation in a sequence of bifurcations of new 
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standing waves from the constant solution (I1.8I1 . The second bifurcation is induced due to the 
positive eigenvalue lo\ in cj(—A), which corresponds to the even eigenfunction in /_ U Iq U /+. The 
corresponding result is described in the following lemma. 

Lemma 3.2. For A < — the constant solution m.8\) is no longer the ground state of the 
constrained minimization problem a for Qq > Qq := ^{L + 27r)n^. Moreover, the constant 
solution is a saddle point of under fixed Q(T) = Qq with one negative eigenvalues for 

Qo £ (Qo> Qo*) negative eigenvalues for Qq > Qq*, where Qq* := \{L + 27r)w^. 

Proof. When the perturbation to the constant solution (jl.8j) is decomposed into the real and imagi¬ 
nary parts U and W, the second variation of = E — AQ is defined by two Schrodinger operators 
L_(_ and L_ as follows: 

(3.1) 6^Ha:={L+U,U)l2 + {L_W,W)l2, 
where 

(3.2) L+ = -A-A-6<^^, 

(3.3) L_ = -A-A-2^>^ 

Because <1>^ is bounded on I- U Iq U the operators L_|_ and L_ are also defined on the domain 
P(A) C if {I- U Iq U IjQ. If is the constant solution (ll.8p . then 

L_ = —A and L_|_ = —A -I- 2A. 

Therefore, the smallest eigenvalue of L_ is located at zero and it is simple. It corresponds to the 
phase rotation of the standing wave $. 

Let us now recall the theory of constrained minimization from Shatah-Strauss m and Weinstein 
|28] . If L+ has only one negative eigenvalue, then is not a minimizer of energy H\. Nevertheless, 
it is a constrained minimizer of energy Ha if a certain slope condition is satisfied. Associated with 
the constraint (5('I') = Qq, we can introduce the constrained if space by 

(3.4) L2:={C/gL2(/„U/oU/+): (C/, 4>)i2 = O} . 

The number of negative eigenvalues of L+ is reduced by one in the constrained space (13.411 if and 
only if ^g(4>) < 0 [27l|28]. Moreover, if ^g(<I>) < 0, the constrained minimizer is non-degenerate 
with respect to the real part of the perturbation U. 

The smallest eigenvalue of L_|_ is negative and the second eigenvalue of L+ is located at -|- 2A. 
Since ^Q(4>) = — (L -|- 27r) < 0, the constant standing wave (|1.8p is a constrained minimizer of 
E{^) if A G (—but it is no longer the ground state if A G (—oo, — 

By Remark 12.3p the eigenvalues flf and ujf in (t(—A) are simple. As a result, operator L_|_ has 
exactly two negative eigenvalues for —^cvi < A < and exactly three negative eigenvalues 

for A < The number of negative eigenvalues of L+ is reduced by one in the constrained 

space Therefore, the constant solution (|I.8p is a saddle point of il'(T) under fixed QPT) = Qq 
with exactly one negative eigenvalue for Qq G {QqiQIq) and exactly two negative eigenvalues for 
Qo ^ Qo*! where Q*q := i(L 2ti)Vl\ and Q*q* := \{L + 27r)a;f. □ 

Remark 3.3. Lemmas 13.II and 13.21 prove the first assertion of Theorem ll.il 

We shall now describe the first (symmetry-breaking) bifurcation of the constant standing wave 
()1.8p at A = or Qq = Qq. We will show that it is a pitchfork bifurcation, which leads to a 

family of positive asymmetric standing waves of the stationary NTS equation (jl.Sp for A < — 

The asymmetric states are the ground states of the constrained minimization problem (II.7p for 
Qo > Qq. The corresponding result is described in the following lemma. 
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Lemma 3.4. Let Aq := — There exists A* G (—oo, Aq) such that the stationary NTS equation 
with A G (A*,Ao) admits a positive asymmetric standing wave which converges to the 
constant solution in the H'^-norm as A —>• Aq. Moreover, there exists Qq G (QqjOo) such that 
the positive asymmetric standing wave is a ground state of the constrained minimization problem 
n for Qo G {Q*o,Q*o). 

Proof. By Proposition 12.11 the eigenfunction U corresponding to the eigenvalue 12^ in a{—A) is 
odd with respect to the central point x = 0, whereas the constant solution (|1.8h is even. Therefore, 
we have the symmetry-breaking bifurcation, which is similar to the one studied in [16]. In order 
to unfold the bifurcation, we study how the odd mode U can be continued as A is defined near 
the bifurcation value Aq := We use the explicit expression for U given by (I2.14[) and the 

characterization of the values of given by Proposition 12.41 By Remark 12.31 is a simple 
eigenvalue in cj(—A). 

Using a simplihed version of the Lyapunov-Schmidt reduction method m, we consider a regular 
perturbation expansion for solutions of the stationary NLS equation (|1.5p near the constant solution 
m cit A. — Aq. Thus, WG Gxpcinci 

(3.5) A = + 0^11 -|- O(a^), 4’(x) = 2^1 a^4’2(x) -|- Ofj2 {p^), 

where a is a small parameter for the amplitude of the critical odd eigenfunction U of the operator 
—A in 'D(A), whereas the corrections O and {4>„}„>2 G T’(A) are defined uniquely under the 
constraints {U,^k)L^ = 0, k > 2. Note that the decomposition (13.51) already incorporates the near¬ 
identity transformation that removes quadratic terms in a and ultimately leads to the normal-form 
equation, derived in a similar context in m- 

At a^, we obtain the inhomogeneous linear equation 

(3.6) (-A-f2?)4>2 = ^Oi(b? + 6C/2). 

Since U is odd, the right-hand-side is even. Thus, the Fredholm solvability condition is satisfied 
and there exists a unique even solution for <I' 2 , which can be represented in the form 

(3.7) 4>2(x) =-^-h 3ni$2(a:), 

where 4*2 is uniquely defined from the linear inhomogeneous equation 

(3.8) {-A-nl)^2 = u^- 
At a^, we obtain another inhomogeneous linear equation 

(3.9) (-A-0?)4'3 = fit/+ 6f2iU4>2+ 21/1 

The right-hand-side is now odd and the Fredholm solvability condition produces a nontrivial equa¬ 
tion for fi: 

(3.10) ^^11^711^2 + 6fii(t/2, 4>2)i2 + 2||t/||i4 = 0. 

Substituting (|3.7I1 into ()3.10p . we obtain 

(3.11) fi||t/||i2 =9fi2([/2,l>2)L2 + ||[/||i4. 

We need to show that the right-hand-side of equation (13.lip is negative, which yields fi < 0. In 
view of the decomposition (13.5p . for a sufficiently small, the new solution represents a positive 
asymmetric standing wave satisfying the stationary NLS equation (11.51) with A < Aq. This would 
imply the first assertion of the lemma. 
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Using the explicit representations (12.141) and (I2.15p . we obtain an explicit solution of the linear 
inhomogeneous equation (13.8p : 


l'2(x) = 


^cos(nix) — [cos(2nix) + 3] , 

5cos(Oi(x - L - tt)) + [cos(2fii(x - L - vr)) - 3] , 


X E Jo, 
X E /+, 


where A and B are constants of integration to be defined, the symmetry <h 2 (~x) = 4>2(x) can be 
used, and the homogeneous sinusoidal solutions are thrown away since sin(7rni) ^ 0. 

Using Kirchhoff boundary conditions (|1.2p . we uniquely determine constants A and B from the 
following linear system of algebraic equations: 


(3.12) 


cos (Till) 

— cos(7rni) 

'A 

sin(Lni) 

2 sin(7rni) 

B 


1 


3 + cos(2Lni ) + cS\tnl] (cos(27rlli) - 3) 
2sin(2Lni) - sin(27rni) 


Using the transcendental equation (|2.15p . we can see that the second entry in the right-hand side 
of (I3.12P is zero. As a result, the second equation of the system (I3.12p yields 


(3.13) 


sin(Lfli)A -|- 2 sin(7rlli)i? = 0. 


Let us prove that the system ()3.12p admits the unique solution in the following explicit form: 


(3.14) 


A = 


■ cos^(Lni 


B = - 


sin^ (LQi ) cos^ (Lfli ) 


’ 2n^cos(7rni) 

Indeed, the constraint (|3.13p is satisfied with the solution ()3.14p . Furthermore, the first equation of 
the system (13.121) is satisfied with the exact solution (j3.14p if and only if the following transcendental 
equation is met: 


(3.15) 


^ sin^(Lni) 1 

cos2(7rlli) 3 


cos(2Lni) + 


sin^(Lni) 


cos(27rni 


= cos^(Lni). 


cos2(7rlli) 

Using the transcendental equation (j2.15p . we rewrite equation (I3.15h in the equivalent form 


(3.16) 


sin2(Lfli) 3 2^ro ^ 

1 - TTTr^ = 

cos^(7riZi) 4 


which is satisfied identically, thanks again to the transcendental equation (j2.15p . 

Using (j3.14p in the expression for $ 2 ; we rewrite the expression for + 911^62 explicitly as 

{ I cos^(Lni) cos(nix) — 2 cos(2nix) — 4, x G Iq, 

cos^( 7 rDi) [“I cos^(Lfli) cos(7rOi) cos(fli(x — L — tt)) + 2 cos(211i(x — L — tt)) — 4] , X G /+. 

After computations of the integrals and simplifications with the help of the transcendental equation 
()2.15p . the right-hand side of equation (|3.1ip is simplified to the form 


(3.17) 


= -3L 



sin(2LQi) 
2 Lni J 


9 nl{u^^ 2 )L^ + \\u\\l, 


— Gtt 


sin^(Lni) 

cos‘^(7rni) 


3sin^(L£li) cos(Lni) 
cos^(7rr2i) 


[l -|- 2 cos^(LQi)] . 


By Proposition 12.41 we have 2LQi < tt for every L, so that every term in (|3.17p is negative. 
Therefore, fl < 0 in (|3.11l) . 

Thus, the first assertion of the lemma is proved. In order to prove the second assertion of the 
lemma, which states that the positive asymmetric standing wave given by the decomposition (13.5p 
is a minimizer of the constrained minimization problem (|1.7I) for we need to compute the 

negative eigenvalues of the operators L_|_ and L_ in (13.2p and (13.31) for <I>. Since L_<I> = 0 and <I> 



























































12 


J. L. MARZUOLA AND D.E. PELINOVSKY 


is positive, the zero eigenvalue of L_ is the smallest eigenvalue of L_. The smallest eigenvalue is 
simple. It corresponds to the phase rotation of the standing wave 

Since <I> bifurcates from the constant solution the operator L_|_ has a unique negative 

eigenvalue and a simple zero eigenvalue at A = Aq. We shall now construct a regular perturbation 
expansion for small |A — Aq] in order to prove that the zero eigenvalue becomes a small positive 
eigenvalue for the bifurcating solution (|3.5p . Therefore, we expand 

L+ = -A-Ql- Qa^hU - - a^Q + OL-(a^). 

The zero eigenvalue of —A — corresponds again to the eigenfunction U. For a sufficiently small, 
we expand the eigenvalue A and the eigenfunction u of the operator L_|_: 

(3.18) \ = a'^A + 0{a^), u{x) = U{x) + aUi{x) + a‘^U 2 {x) + 


where the corrections A and {Un}n>i £ are defined uniquely under the constraints {U, Uk)L 2 = 

0 , A; > 1 . At a, we obtain the inhomogeneous linear equation 

(3.19) (-A - nl)Ui = 

which has the unique even solution Ui = 611 1 ^ 2 - At o^, we obtain the inhomogeneous linear 
equation 

(3.20) (-A- n‘l)U2 = miUUi + 6ni^2U + -m + au. 

With the account of (|3.7p . (|3.1ip . and Ui = 6 ni$ 2 , the Fredholm solvability condition yields 


(3.21) A\\U\\l, = -iml{U\^2)L^-2\\U\\\,. 

Comparison with (13.lip yields A = —211. Since we have already proved that H < 0, we obtain A > 0, 
so that for |A — Aq] sufficiently small, the operator L_|_ has a small positive eigenvalue bifurcating 
from the zero eigenvalue as A —>■ Aq. Thus, the operator L+ has only one simple negative eigenvalue 
for A < Aq. 

It remains to show that Q computed at the positive asymmetric standing wave $ given by the 
decomposition ()3.5I) is an increasing function of the amplitude parameter a. In this case, the slope 
condition ^Q(‘h) < 0 holds and the only negative eigenvalue of operator L+ is removed by the 
constraint in defined by (13.411 . From (|1.6I1 . (13.51) . (13.61) . and (13.81) . we obtain 

Q(<h) = g5 + a2(lli(l,$2)L2 + ||[/||2,)+0(a3) 

= Ql-a^ {a{L + 2Ti) + 2 \\U\\l 2 )+0{a^). 


We observe 11(L + 27r) + 2||[/|||2 < 0 if and only if 


(3.22) 

Using (13.17P and 


(L + 27r)(911?([/2, 3.2)i2 + \\U\\\,) + 2\\U\\\2 < 0. 


sin^(LHi 


^2 = L — sin(Llli) cos(Llli) + 2- ^ [tt + sin( 7 rHi) cos( 7 rHi)] 

COS I TT' u ) 


= L + 27r 


sin^(LHi) 


cos^(7rHi) ’ 

the left-hand side of (|3.22p can be written as 

-3(L + 2 ^) (l (1 - [1 + 2cos2(M.)] } 


+2 


L + 27r 


sin^(Lni) 

cos2(7rr2i) 
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We regroup these terms as the sum / + // + III, where 

,4 


I = -3(L + 27r) 


3L sin^(LOi) 

4 cos^(7rr2i) 


+ 2 


L + 27r 


sin^(Lr 2 i) 

cos2(7rOi) 


II = --{L + 2tt)L 


1 - 


sin( 2 Lr 2 


i) 


III = -3{L + 27r) 


2 LOi 

sin(Lr2i) cos(Lrii) f sin^(Lri 


r 2 i 


008^(71^2 


After expanding the brackets, the first term becomes 


1 


1 


I = --L^ - -ttL 


sin^(Lr2i) sin‘^(-Lr2i 

9-16-+ 12 ^ 


^ [l + 2cos2(Lfii)] - . 

2 sin'^(L 12 i) 


cos2(7rr2i 


008*^(71121 


— 471^ 


008 *^( 71121 ) ’ 


where every term is negative because f{x) = 9 — 16x^ + 12x^ > ^ > 0. Hence, I < 0. Furthermore, 
// < 0 holds. Since L12i < ^ by Proposition 12.41 then III < 0 if and only if 5 ( 12 i) > 0, where 

Using equation (I3.16p . we rewrite the left-hand-side as follows: 


5 (^ 1 ) 


1 — - cos^(L12i) 


[1 -|- 2cos^(L12i)] — 


- [1 — cos^(L12i)] [1 -|- 6 cos^(L12i)] 


3 

4 


from which it follows that 5 ( 12 i) > 0 for every L > 0. Thus, III < 0, so that Q(4>) is an increasing 
function of the amplitude parameter a. 

Since ^Q(^) < 0, the operator L+ does not have a negative eigenvalue in the constrained space 
so that 4> is a ground state of the constrained minimization problem (11.71) . The statement of 
the lemma is proved. □ 


Remark 3.5. Asymptotic expansions similar to the ones used in the proof of Lemma 13.41 can be 
developed for the second (symmetry-preserving) bifurcation of the constant standing wave (II.Sp at 
A = —\oj\ or Qq = Qq* := ^(L-|-27r)a;f, where the positive eigenvalue (jj\ in f 7 (—A) corresponds to 
the even eigenfunction. As a result of this bifurcation, a new family of positive symmetric standing 
waves exists for A < —\oJi or Qq > Qq*. The positive symmetric wave is not, however, the ground 
state of the constrained minimization problem (11.71) for Qq > Qq*, because the operator L+ has 
two negative eigenvalue and a simple zero eigenvalue at A = — Even if the zero eigenvalue 
becomes small positive eigenvalue for A < and if one negative eigenvalue is removed by a 

constraint in L^, there operator L_|_ still has one negative eigenvalue in L^. 

Remark 3.6. Lemma[33]and the result stated in Remark l3.5l prove the second assertion of Theorem 

o 


4. Proof of Theorem O 

Here we study standing wave solutions of the stationary NLS equation (11.51) in the limit A —— 00 . 
The standing wave solutions are represented asymptotically by a solitary wave of the stationary 
NLS equation on the infinite line. We use the scaling transformation 

(4.1) 4>(x) = |A| 2 'I'( 2 ;), 2 ; = |A| 2 x, 
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and consider the positive solutions of the stationary NLS equation (|1.5I1 for A < 0. The stationary 
problem can then be written in the equivalent form 

(4.2) - + = 0, zEJ_UJoUJ+, 

where is the Laplacian operator in variable z and the intervals on the real line are now given by 

J- :=[-{L + 2TT)fj,,-Lfj], Jo := J+:= [L/r, (L + 27r)/r] , 

with /i := |A| 2 . The Kirchhoff boundary conditions (II.ip and (II.2|] are to be used at the two 
junction points. 

The stationary NLS equation — A^T + T — 2T^ = 0 on the infinite line is satisfied by the solitary 
wave 

(4.3) 'I'oo(-2) = sech(z), z G R. 

To yield a suitable approximation of the stationary equation (14.2p on the dumbbell graph J_ U 
Jo U J+, we have to satisfy the Kirchhoff boundary conditions at the two junction points. Two 
particular configurations involving a single solitary wave will be considered below: one where the 
solitary wave is located in the central line segment and the other one where the solitary wave is 
located in one of the two loops. As follows from numerical results reported on Figures [H [71 El and [9] 
below, these configurations are continuations of the two families of positive non-constant standing 
waves in Lemma 13.41 and Remark 13.51 

4.1. Symmetric solitary wave. We are looking for the symmetric standing wave 

(4.4) ^'o(-^) = ^o(^), 2 € Jo, ^'-(-2) = ^+(z), ^ G J+. 

We will first provide an approximation of the solitary wave with the required Kirchhoff boundary 
conditions by using the limiting solitary wave (|4.3p . Then, we will develop analysis based on the 
fixed-point iterations to control the correction terms to this approximation. The following lemma 
summarizes the corresponding result. 

Lemma 4.1. There exist /tq > 0 sufficiently large and a positive fi-independent constant C such 
that the stationary NLS equation J[ ) for p, G (/io,oo) admits a symmetric standing wave T near 
Too satisfying the estimate 

(4-5) l|4'- 4 'oo||l°o(j_uJouj+) < 

Proof. The proof consists of two main steps. 


Step 1: Approximation. We denote the approximation of T in Jq by Go and set it to 


(4.6) 


Go{z) = sech(z), z ^ Jq. 


The approximation of T in J_|_, denoted by G+, cannot be dehned from a solution of the linear 
equation G'f — G+ = 0 because the second-order differential equation does not provide three 
parameters to satisfy the three Kirchhoff boundary conditions: 

(4.7) G+{Lp) = G+{{L + 27 t)p) = Go(L/r), GV(L^) - GV((L + 27r)/r) = G'(L/r). 


Instead, we construct a polynomial approximation to the Kirchhoff boundary conditions, which 
does not solve any differential equation. Using quadratic polynomials, we satisfy the boundary 
conditions (j4.7p for the solitary wave (14.6p with the following approximation 


(4.8) 


G+{z) = sech(L^) 


1 + 


tanh(L/r) 

47r/i 


{z — Lp.){z — {L + 27r)/i) 
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Note that the maximum of G+ occurs at the middle point of J+ at z = (L + tt)//, and for sufficiently 
large we have 

(4.9) 

for a positive /i-independent constant C. Also, the first and second derivatives of G+ do not exceed 
the npper bound in (14.91) . 

Step 2: Fixed-point arguments. Next, we consider the correction terms to the approximation 
G in (j4.6p and (14.8p . Using the decomposition ih = G + ip, we obtain the persistence problem in 
the form 

(4.10) L^iP = Rqs{G) + N{G,iP), 

where := —+ 1 — 6G{z)^, Res(G) := A^G — G + 2G^, and N{G,ip) := QGip"^ + 2ip^. Since 
the approximation G satisfies Kirchhoff boundary conditions, which are linear and homogeneons, 
the correction term ip is required to satisfy the same Kirchhoff boundary conditions. The residnal 
term Res(G) is snpported in J_ and J+ and it satisfies the same estimate as in (14.9h . Transferring 
this estimate to the norm, since the length of J+ grows linearly in /x, we have for all snfficiently 
large /i, 

(4.11) l|Res(G)|U2(j_ujouJ+) < 

where the positive constant G is /U-independent. 

The operator is defined on L?{J- U Jq U J+) with the domain in T){Az), that incorporates 
homogeneous Kirchhoff boundary conditions. As /x —oo, the operator converges pointwise to 
the operator 

(f 

L^:= -— + 1- 6sech2(z) : ^'^(M) ^ 

which has a one-dimensional kernel spanned by whereas the rest of the spectrum of Loo 

includes an isolated eigenvalue at —3 and the continuous spectrum for [l,oo). Since the operator 
Loo is invertible in the space of even functions, it follows that the operator L^ : L>(A^) —)• L^(J_ U 
Jo U J+), is also invertible with a bounded inverse on the space of even functions if /U is sufficiently 
large. In other words, there is a positive /x-independent constant G such that for every even 
/ G L^{J- U Jo U J+) and sufficiently large /x, the even function L~^f satisfies the estimate 

(4.12) \\I^Pi^f\\H^{j.uJoUJ+) < C'II/IIl2(j_uJouj+)- 
Hence we can analyze the fixed-point problem 

(4.13) ip = L;1 [Res(G) + N{G, iP)] , iP e V{A,) 

with the contraction mapping method. By using (14.111) . (I4.12p . and the Banach algebra properties 
of U Jo U J+) in the estimates of the nonlinear term N{G,u), we deduce the existence of a 

small unique solution ip G T>{Az) of the fixed-point problem (14.131) satisfying the estimate 

(4.14) II'*/’IIh2(j_uJouj+) < 

for snfficiently large /x and a positive /U-independent constant G. By the constrnction above, there 
exists a solntion T = G -|- V' of the stationary NTS equation (|4.2p that is close to the solitary wave 
(14.3p placed symmetrically in the central line segment. The estimate (j4.5p is obtained from (|4.14|) 
by Sobolev’s embedding of H^{J- U Jo U J+) to L°°{J- U Jo U J+). □ 

Remark 4.2. The method in the proof of Lemma l4.ll cannot be nsed to argue that T is positive, 
although positivity of T is strongly expected. In particnlar, G+ is not positive on J+ and the 
correction term x/;_|_ is comparable with G+ in J+. Similarly, the approximation G of the solntion T 
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is not unique, although for every G, there exists a unique correction if; by the contraction mapping 
method used in analysis of the fixed-point problem (|4.13p . We will obtain a better result in Lemma 
K7\ with a more sophisticated analytical technique in order to remove these limitations of Lemma 

KT\ 

4.2. Solitary wave in the ring. We are now looking for an approximation of the solution Ik of 
the stationary NLS equation (14.2p . which represents as —>■ oo a solitary wave residing in one of 
the rings, e.g. in J+. Because of the Kirchhoff boundary conditions in J+, the approximation of ik 
in J_|_, denoted by G+, must be symmetric with respect to the middle point in J_|_. Therefore, we 
could take 

(4.15) G+{z) = sech(z — L/j. — TTfi), z £ J+. 

However, the method used in the proof of Lemma 14.11 fails to continue the approximation (|4.15p 
with respect to finite values of parameter fi. Indeed, the linearization operator : 'D(Az) —)• 

L^(J_ U Jo U J+) defined on the approximation G+ has zero eigenvalue in the limit —)• oo, which 
becomes an exponentially small eigenvalue for large values of /r. Since no spatial symmetry can 
be used for the correction term to G+ on J_|_ because of the Kirchhoff boundary conditions in 
J+, it becomes very hard to control the projection of 4^+ to the subspace of V(Az) related to the 
smallest eigenvalue of L^. 

To avoid the aforementioned difficulty and to prove persistence of the approximation ()4.15p . we 
develop here an alternative analytical technique. We solve the existence problem on J_|_ in terms 
of Jacobi elliptic functions with an unknown parameter and then transform the existence problem 
on J_ U Jo with the unknown parameter to the fixed-point problem. After a unique solution is 
obtained, we define a unique value of the parameter used in the Jacobi elliptic functions. The 
following lemma summarizes the corresponding result. 

Lemma 4.3. There exist fJ-o > 0 sufficiently large and a positive ^-independent constant C such that 
the stationary NLS equation for fi € {ixo,oo) admits a unique positive asymmetric standing 
wave 'k given by 

(4.16) 

and satisfying the estimate 

(4-17) ||4'o||h2(Jo) + I|4'-II//2(j_) < Ce~'^^, 

where dn(^; k) is the Jacobi elliptic function defined for the elliptic modulus parameter k G (0,1). 
The unique value of k satisfies the asymptotic expansion 

(4.18) [1 + as ^ ^ oo. 

Remark 4.4. It is well-known (see, e.g.. Lemma 2 in [22] for similar estimates) that if /r —)• (X) and 
k ^ 1 according to the asymptotic expansion (I4.18P , then the dnoidal wave (I4.16P is approximated 
by the solitary wave (|4.15l) on J+ such that 

(4.19) ||'k+— < Ge 

where the positive constant G is //-independent. Therefore, for sufficiently large /t, we have justified 
the bound 

ll'k - 'koo(- - Lp.- p)\\l°°{G-vjGoVJG+) < Ce~'"^, 
for the asymmetric standing wave of Lemma 14.31 
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Proof. The proof of Lemma 14.31 consists of three main steps. 

Step 1: Dnoidal wave solution. The second-order differential equation (14.211 is integrable 
and all solutions can be studied on the phase plane The trajectories on the phase plane 

correspond to the level set of the first-order invariant 

2 

_ \J/2 -(- \J/4 _ 

The level set of I is shown on Figure [3l There are two families of periodic solutions. One family 
is sign-indefinite and the corresponding trajectories on the phase plane surround the three 

equilibrium points. This family is expressed in terms of the Jacobi cnoidal function. The other 
family of periodic solutions is strictly positive and the corresponding trajectories on the phase plane 
(T,'!'') are located inside the positive homoclinic orbit. This other family is expressed in terms of 
the Jacobi dnoidal function. 

0.8 r 


0.4 H 


© 0 


-0.4 


_0,8l-!-1-!-^^-1 

-1.5 -1 -0.5 0 0.5 1 1.5 

<I> 




(4.20) 


I := 


dz ) 


Figure 3. Level set (I4.20jl on the phase plane (T, T'). 

Because of the Kirchhoff boundary condition T_|_(L/i) = ^^{Lfi + 2TTfi), we consider a trajectory 
on the phase plane (T, T'), which is symmetric about the middle point in J_|_ at z = (L-|-vr)/i. Since 
the trajectory is supposed to converge to G+ given by (j4.15jl as /r —>■ oo, we select an incomplete orbit 
with 'LY(L/r) = + 27r/i). For the trajectory inside the homoclinic orbit, the corresponding 

solution of the first-order invariant (j4.20p is given by the exact expression (j4.16p . We now define 

= dn ( ^ , 

+ ^ V2^ ) 

f'UM = cn 

and consider the range of the values of k for which q{k, fi) > 0. Since cn(^; k) vanishes at .^ = K{k), 
where K{k) is the complete elliptic integrals of the first kind, we can define from the unique 

root of the equation 

(4.23) TTfr = ^2-k^K{k^). 

Therefore, q{k^{fi),fi) = 0. On the other hand, we have 

(4.24) qik,fi) —)• q*{fJ-) ■= tanh(7r/i)sech(7r^) = 2e“^^ -|- as /c —)• 1. 


(4.21) p{k,n) := 

(4.22) q{k,ij) : = 
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Since (see 8.113 in [12]) 

(4.25) K{k) = log +0 {{I-k^)\\og{l-k^)\) as k ^ I, 

the root A:* (/r) of the transcendental equation (|4.23D satisfies the asymptotic expansion 

\/l — = 46”^^^ + 0{e~^^^) as /r — )• oo. 

Thus, the entire interval 1) is exponentially small in terms of large /r. With this asymptotic 

in mind, we compute the limiting values of the function p{k,^) for k G {kif{fi), 1). At one end, we 
obtain 

(4.26) p{k^{p),p) = ^ = 4e“^^ + as oo, 

\/2 - k ^ 

whereas at the other end, we obtain 

(4.27) p{k,p) p*{p) ■= sech(7r^) = 2e~^^ + as A —)• 1. 

Figure m shows the dependencies of p and q versus k in {k^{p), 1) for a particular value p = 2. 
The graph illustrates that q{k,p) is a monotonically increasing function with respect to A in /c G 
{k^{p), 1) from 0 to q*{p), whereas p{k,p) is a monotonically decreasing function from p{k^{p),p) 
to p*{p). 



Figure 4. The graphs of p (solid blue), q (solid red), and 2q (dashed red) versus k 
in {k^{p), 1) for p = 2. The dotted line shows a graphical solution of the equation 
p = 2q at ko{p). 


Step 2: Fixed-point arguments. By substituting the explicit solution (I4.16[) to the stationary 
NTS equation (|4.2p . we close the system of two second-order differential equations 


(4.28) 


r + 4'o - 24^^ = 0, zGJo, 

\ -4'" -h 4'_ - 24'?. =0, z G J_, 


with five boundary conditions 

( 4'_(-L/i) = 4'_(-L/r - 27r^) = 4'o(-Tp), 

(4.29) 4''_(-L/r)-4/'_(-L/r-27r^) = 4''(-L/i), 

[ 4'o(T/r) = p{k, p), %{Lp) = 2q{k, p), 

where p{k, p) and q{k, p) are explicit functions of an unknown parameter k defined in the interval 
{k^:{p),l). Thus, although the five boundary conditions over-determine the system of differential 
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equations (|4.28p . the parameter k can be used to complete the set of unknowns. For this step, we 
will neglect the last boundary condition in (14.291) given by = 2q{k,fi). 

Let us define the approximation of the solution to the system (14.281) . denoted by (G_,Go) by 
the solution of the linear system 


f ~Go + Go — 0) z £ Jq 
\ -G'L + G_ = 0, z G J_ 


subject to the four boundary conditions 


r G-{—L^) = G-{—Lfj, — 2Trfi) 
(4.31) < G-i-L^i)-G'_i-Ln-2Trn) 

1 Go{Ln) = p{k,n). 


Go{—Lp,), 

G'o(-L^), 


The boundary-value problem has the unique solution 

COsh(z -|- Lfi + TTfi) 


(4.32) 

(4.33) 


G-{z) =p{k,fj.) 
Go{z) =p{k,p) 


cosh(7r/x) cosh(2L/x) -|- 2 sinh(7r;u) sinh(2L/i) ’ 
cosh(7r/x) cosh(z + Lp) + 2 sinh(7r/i) sinh(z + Lp) 


cosh(7r/i) cosh(2L/i) + 2 sinh(7r/i) sinh(2L^) 

It follows from the explicit solution (14.321) and (14.331) that for p > 0 sufficiently large, we have 
(4.34) \\G-\\hHj.) < Gp{k,p)e-^^>^, WGoWh^Jo) < Gp{k,p), 


where p{k,p) = 0(e ’^^) as p ^ oo for every k G {k^{p),l) and G is a positive //-independent 
constant. 


Remark 4.5. Compared with the estimates (j4.9h and (14.lib , where we are losing p^G between the 
L°° and bounds, the estimate (14.341) works equally well in L°° and thanks to the integration 
of the explicit solutions (|4.32p and ()4.33p . 


Applying the decomposition 


T_(z) = G_(z) -ip-iz), To(z) = Go(z) 'ipo{z), 


we obtain the following persistence problem 

qr'i / —i’o + V'o = 2(Go -I- V’o)^; ^ £ Jo, 

^ ^ \ TV'-= 2(G_+V’-)^ ^gJ_, 

subject to the four homogeneous boundary conditions 

( i’-i-Lp) = '4^-{-Lp - 27r//) = ^poi-Lp), 

i^'_{-Lp) - i^'_{-Lp - 2Trp) = V’o(-L//), 

V’o(^h) = 0- 


Since 1 — is invertible on L?‘{J- U Jo) with a bounded inverse in H^{J- U Jo) due to the four 
symmetric homogeneous boundary conditions (I4.36P and the inhomogeneous term is estimated by 
using the bound (I4.34h . a contraction mapping method applies to the fixed-point problem (I4.35h . 
As a result, there exists a small unique solution for and 'i/jq of the persistence problem (|4.35l) 
satisfying the estimate 

(4-37) IIV'-IIh2(j_) + IIV'oIIh2(Jo) < Gp{k,p)^, 


where p{k,p) = 0(e as // —>■ oo for every k G {k^{p),l) and the positive constant G is p- 
independent. 
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Step 3: Unique value for the parameter k in the interval {k^{pL), 1). It remains to satisfy 
the fifth boundary condition in (j4.29p . which can be written in the following form 

(4.38) 2q{k, fj.) = G'oiLfi) + = p{k, p) + as p^oo, 


where we have used the exact result (14.331) yielding 

m) f 1 + 


2sinh(7r/i) — cosh(7r/i) 


cosh(7r/r) cosh(2L^) + 2 sinh(7r/i) sinh(2L/i) 


as well as the estimate (j4.37l) with the account that p{k,p) = 0(e as /i —>■ oo for every 
k e (fc*(p), 1). 

Figure 0] illustrates graphically that q(k, p) is monotonically increasing with respect to k in the 
interval {k^{p), 1) from 0 to q*{p) = 2e~'"^ + whereas p{k, p) is monotonically decreasing 

with respect to k in the interval {k^{p), 1) iiom p{k^{p), p) = 4e“^^+ to p^{p) = 2e~'"^ + 

0(e-3^^). Therefore, there exists exactly one solution /co(m) £ {k*{p)-, 1) of the equation (|4.38l) . 

In order to prove monotonicity of p and q in k, as well as the asymptotic expansion (I4.18p . we 
study the rate of change of the functions p and q with respect to k. From the explicit expression 
(I4.2ip . we obtain 


(4.39) 


dkpik, p) = 


dn 


^/{2-k‘^f \V2^ 


Tip 


TT 


pk^ 


rsn 


np 


+ 


(2-P)2 

1 


k I cn 




5fcdn (^;/c) 




irp 

V2-k^ 



For every k E {k^,{p), 1) and sufficiently large p, the first term in (I4.39P is exponentially small of 
the order of 0{e~'^^), whereas the second term is larger of the order of 0{pe~'"^). Nevertheless, 
we show that the last term in (j4.39p is dominant as it is exponentially large as /i —)• oo. To show 
this, we use the following result proved in Appendix A. 


Proposition 4.6. For every ^ E M, it is true that 

(4.40) sn(^; 1) = tanh(0, dksn{^; 1) = [sinh(0 cosh(0 - C] sech2(^), 

(4.41) cn(,^; 1) = sech(^), dkcn{^; ~ ^ [sinh(^) cosh(^) — tanh(^)sech(^), 

(4.42) dn(^; 1) = sech(^), 9fcdn(^; 1) = —^ [sinh(^) cosh(^) + tanh(^)sech(^). 

Moreover, if p is sufficiently large, then for every ^ E {0,Trp) and every k E {k^{p),l), there is a 
positive p-independent eonstant C sueh that 

(4.43) |4sn(^; k) - 4sn(C; 1)| + |4cn(C; k) - dkcn{^; 1)| + |4dn(^; k) - 4dn(^; 1)| < Cpe~^>^. 
From (I4.42D and (|4.43ll . we obtain the dominant contribution of (j4.39D for every k E {k^{p), 1): 

(4.44) dkpfk, p) = + 0{pe~'^^) as p ^ oo. 

Similarly, we differentiate (|4.22p in k, use (|4.40p . (j4.41l) . and (I4.43p . and obtain the asymptotic 
expansion for every k E {k^,{p), 1): 

dkq{k, p) = + 0{pe~'"^) as p ^ oo. 


(4.45) 
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It follows from (|4.44l) and (14.451) that p{k, fi) and q{k, fi) are monotonically decreasing and increasing 
functions with respect to /c as /c —>■ 1, in agreement with the behavior on FigurelH Furthermore, the 
algebraic equation (I4.38P can be analyzed in the asymptotic limit of large p. Indeed, multiplying 
()4.38l) by we obtain 

4g-2^/. + i(A; _ 1) + {k - 

(4.46) = -^{k-l)+ {k - g-27rM-4LM)^ 

where remainder terms are all smooth in their variables. By the Implicit Function Theorem, we 
obtain the unique root of the algebraic equation ()4.46p denoted by ko{p). The root satisfies the 
asymptotic expansion 

koip) = 1 - [1 + , 

which justifies the asymptotic expansion (|4.18p . Furthermore, the bound (I4.17P follows from esti¬ 
mates (j4.34p and (j4.37p . 

Finally, because the perturbation term i/^o) is triply exponentially small, whereas the leading- 
order approximation (G_,Go) is exponentially small and positive, we deduce that (4'_,'I'o) is 
positive on J_ U Jq. From the exact representation (|4.16l) . we also know that 'I'_|_ is positive on 
J_i_. Thus, the asymmetric standing wave is positive on J_ U Jq U J+. The proof of the lemma is 
complete. □ 


The same method in the proof of Lemma 14.31 can be applied to construct the symmetric solitary 
wave described in Lemma O However, because of the Kirchhoff boundary conditions, we need 
to take the symmetric orbit outside of the homoclinic orbit on Figure [3l The following lemma 
summarizes the corresponding result. 


Lemma 4.7. There exist /tq > 0 sufficiently large and a positive fi-independent constant C such 
that the stationary NLS equation for p G {po, oo) admits a unique positive symmetric standing 
wave 'k given by 


(4.47) To(z) 

and satisfying the estimate 


\/2A;2 - 1™ (v2A:2 - l’ 


Z G Jq 


(4.48) 


^+\\h2{j+) + II^-IIh2(J-) - 


where cn(|; A:) is the Jacobi elliptic function defined for the elliptic modulus parameter k G (0,1). 
The unique value for k satisfies the asymptotic expansion 


(4.49) 


■\/l — kJ 


-^e ^^\l + 0{pe '^^^)] as p ^ oo. 


Proof. We only outline the minor differences in the computations compared to the proof given in 
Lemma 14.31 For the trajectory outside the homoclinic orbit, the corresponding solution of the 
first-order invariant (j4.2np is given by the exact expression (j4.47l) . We now define 


(4.50) 

(4.51) 


p{k,p) 

q{k,p) 


'^oiLp) = 


k 


icn 


( Lp 


V2A:2 - 1 \^J2k‘^ - 1 


]k 


-<{Lh) = 


k 


-sn 


Lp 


2k‘^-l V\/2A:2 - 1 


k dn 


Lp 


V2k‘^ - 1 
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The trajectory is already even in z. We consider the range of the values of k for which p(k, /x) > 0. 
Therefore, k is defined in {k^{fi), 1), where /c*(/x) is the root of the algebraic equation 

/xL = y^2kl^K{k,), 

which is expanded asymptotically as 

\/l — as /U —)• oo. 

Again, the interval (A:*(/x), 1) is exponentially small as ^ oo. 

Figure [5] shows the dependencies of p and q versus k in {k^(p), 1) for a particular value p = 2. 
The graph illustrates that p{k, p) is a monotonically increasing function with respect to k from 0 
at fe = k^{p) to 

p*{p) := sech(L;u) = 2e~^^ + 

as /c ^ 1, whereas q{k,p) is a monotonically decreasing function in k from 

q{K{p),p) = =4.e-^^ + 0{e-^^n 

sX k = k^,{p) to 

q^{p) := tanh(L|u)sech(L|u) = 2e~^^ + 

as /c —1. Again, there is a unique root k = ko{p) in {k^{p), 1) of the algebraic equation 

q{k, p) = 2p{k, p) + as p^oo. 

The bound (I4.48p . the asymptotic expansion ()4.49p . and the positivity of the symmetric wave are 
proved by similar estimates to those in Lemma 14.31 □ 



Figure 5. The graphs of p (solid blue), q (solid red), and 2p (dashed blue) versus 
k in {k^{p), 1) for ^ = 2 and L = vr. The dotted line shows a graphical solution of 
the equation q = 2p aX ko{p). 


Remark 4.8. Lemma [4.71 gives an improvement of Lemma [4.11 since the symmetric standing wave 
T is now proved to be strictly positive and the estimate (|4.5I) is now improved to be 

(4-52) 11'^'“ ^oo||l°°(j_uJoUJ+) ^ C'e 

where the positive constant C is /x-independent. 
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4.3. The symmetric standing wave as a ground state. There exists a simple argument why 
the ground state of the constrained minimization problem (11.71) is represented by a single solitary 
wave in the asymptotic limit A ——oo. Indeed, computing asymptotically Eq and Qq at the 
representation (jl.lOl) . we obtain 


(4.53) 


2 3 1 

( 5 o~ 2 |A| 2 , as A —)•-oo. 

O 


For N solitary waves packed in the same graph /_ U Iq U /+ at different points, we estimate Qq and 
Eq roughly by multiplying (j4.53|) by N, so that if Qq is preserved, then 


(4.54) 


|A| 2 ~ 


Qo 

2N 


En ~ — 


Q\ 


i2Ar2- 

Therefore, the standing wave of minimal energy as A ^ —oo corresponds to A^ = 1. Two single 
solitary waves are given by Lemmas 14.31 and 14.71 The following lemma specihes that the symmetric 
solitary wave of Lemma 14.71 is a ground state of the constrained minimization problem (11.71) . 

Lemma 4.9. For sufficiently large negative A, the symmetric wave of Lemma has a smaller 
Qo o.t a fixed A compared to the asymmetric wave of Lemma \4-3\ Consequently, it has smaller Eq 
at fixed Qo for sufficiently large Qo ■ 

Proof. By the scaling transformation (14.ip . we have Qq = p,Qo, where 

Qo = ll^lli2(j_ujouJ+) and = |A|^/2_ 

We represent Qo at the asymmetric wave of Lemma 14.31 as the sum of three terms 

Qo = ll^+lli2(J^) + Il^0|li2(jg) + Il^-Ili2(j_)- 

By the estimate ()4.17p . there exists a positive //-independent constant C such that for sufficiently 
large //, we have 

(4.55) lltoll^tj.) + I|4 '-II!»(j_) < Ce-^'e 

On the other hand, by the explicit expression (I4.16|) . we have 


I 


Ii2(j+) 


i 


7rfij\/2—k‘^ 


(4.56) 


V2-A;2 

•\/2 — J(^ 


dn{Q, k)‘^df, 


rK{k} 


dn(^; k)‘^df + 


Tr^/^/2—k'^ 


/•TTfi/y 

JKik) 


dn(^; k.)‘^df^. 


V2 - k^ 

It follows from the asymptotic expansions (|4.18l) and (14.2511 that the second term has the following 
asymptotic behavior as // —)> oo: 


(4.57) 


'Kfl/y/2—k'^ 


JKik) 


dn(^; k)‘^df 


/*71 

J 77 L 


dn{Q, kfdf = 0{e-^^>^), 

/ 7r/4+C?(/xe“'^^) 

since dn(^; kfi = 0{e~‘^'"^^) for every ^ = 0{p) and k G {k^{p), 1) as /t —>■ oo. Therefore, the second 
term in ()4.56l) is comparable with the other remainder terms in (|4.55l) . We will show that the first 
term in (I4.55P has the larger value as // —)• oo. We recall that (see 8.114 in |12] 1 

rK{k) 

dn(,f; k)‘^df 


f 

E{k) := / 
Jo 


1 


= 1 + 1(1-^) 


log 


VI - k'^ 


-|-O ((1 — A;^)^| log(l —/c^)!) as A; —)• 1, 
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where E(k) is a complete elliptic integral of the second kind. As a result, we obtain 


V2 - 


l-Kik) .. 

dn{^;= 2 + 


vr/i + log(Vs) - - 




Combining this estimate with (14.551) . (I4.56p . and (I4.57p . we obtain for the asymmetric wave of 
Lemma 14.31 that 

(4.58) = 2|A|^/2 ^ ^^|^|g-2.|A|V2 ^ ^ |'|^|l/2g-2^|A|V2^ 


as A —)• —oo. 


We now report a similar computation for the symmetric wave of Lemma 14.71 We represent Qq 
for the symmetric wave as the sum of two terms 

Qo = ll^o|li2(jg) + 2||^'+||^2 (j_^). 

By the estimate p4.48p . there exists a positive /r-independent constant C such that for sufficiently 
large /i, we have 


(4.59) 


1^' 


+ IIL2(J+) - 




On the other hand, by the explicit expression p4.47l) . we have 
2^2 —1 


(4.60) 


V2P^ 

2 

V2A:2 - 1 


li 


K(k) 


, /> , x2 jA 2(1 - k‘^)^L 

dn(|; hydi - ^ ’ 


2/c2 - 1 


dn(^; A:)2d^ — 


2(1 — k‘^)^L 

2k^-l ' V2A:2 - 1 


+ 


^ pbL j yj‘Ik^— \ 

fK(k) 


(hi{^]k)‘^d^. 


By the same estimate as in p4.57D . the last term in (|4.6Up is comparable with the estimate (|4.59p . 
whereas the other two terms give a larger contribution. We now compute these terms explicitly 


2 16 „„ 
2e-i 3' 


L/i + log(V3) - - 


+ C>(/ie 


—4L/i —2L/i—477/1 


)■ 


Combining this estimate with (|4.59l) and (I4.60p , we obtain for the symmetric wave of Lemma 14.71 
that 


(4.61) 


gr = 2|A| 


i/2_^i|A|e-2^AI^'^+o(|A| 


l/2^-2L|A|i/2\ 


as 


A^ - 


oo. 


For sufficiently large |A|, we have gg^™ < 2 |AV 2 < which proves the first assertion of the 

lemma. 

We shall now prove that this estimate for gg at a fixed A can be transferred to the similar 
estimate for Eq at a hxed gg. This is done from the variational principle for the standing wave 
solutions of the stationary NLS equation (II.5p : 


which implies that 
(4.62) 


dEf) _ dQ g 

dEo 


dQo 


= A, 
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where A is supposed to be expressed from Qq. Since Qo = 0(|A|^/^) as A —>■ —oo, the dependence 
A e-)• (5o is a decreasing diffeomorphism, which can be inverted. Indeed, for the symmetric wave of 
Lemma 14.71 we have 

|A|^/^ = 2^0 + 2 AQqC + 0{Qoe g^g 

so that (I4.62P implies that 


1 rQo 


+ 0(Qoe-^^°) 


dQo <0 as Qo ^ oo. 


On the other hand, for the asymmetric wave of Lemma 14.31 we have 
|A|V2 = 

so that (j4.62p implies that 


|A|^/2 = ^Qo - ^vrQge-Qo + oiQoe-^^^), as Qo ^ oo, 


pasym ^ 




Qo 

+ 0O _ 


-^T^Qle-^^^ + 0{Qoe-^^°) 


dQo >0 as Qo —)• oo. 


Therefore, for sufficiently large Qo, we have < —^Qq < which proves the second 

assertion of the lemma. □ 


Lemma 4.10. There exists Q” G (Qq*,oo) such that the symmetric wave is the ground state of 
the constrained minimization problem ra for Qo G (Qq*,oo). In particular, the second eigenvalue 
of the linearization operator L+ at the symmetric wave of Lemma E3 is strictly positive. 

Proof. By using the scaling transformation (14.111 . we transform the linearized operators L_|_ and L_ 
given by (|3.2p and (13.3p to the form L± = |A|£±, where 

£+ = -A, + 1 - 64/2, 

£_ = -A^ + 1 - 24/2^ 

where both operators are defined on the domain P(A^) in iffj- U J_ U J+). We consider the 
symmetric standing wave 4/ given by Lemma 14.71 

Since £_4/ = 0 and 4/(z) >0 for every z G J_U J_U J+ by Lemma HTTl the operator is positive 
dehnite. Therefore, we only need to show that the operator has a simple negative eigenvalue and 
no zero eigenvalue. It is clear that is not positive definite because (£+4/, 4/)j;^2 = —4||4/||^4 < 0. 
In the limit /r —oo, converges pointwise to the operator 

f]2 

+ 1 - 6sech2(z) : H^{R) 

which admits a simple negative eigenvalue and a simple zero eigenvalue. Therefore, we only need 
to show that the simple zero eigenvalue of becomes a positive eigenvalue of for large but 
finite p.. 

We note that £_|_4/'(2:) = 0, although 4/'(z) does not satisfy the Kirchhoff boundary conditions 
in P^Az). To correct the boundary conditions, we write an eigenfunction U G P{Az) of the 
eigenvalue problem i^U = XU in the product form U{z) = a{z)^'{z). The amplitude function 
a : J_ U J_ U J+ —)• M ensures that the eigenfunction U satisfies the Kirchhoff boundary conditions. 
If 4/ is even, then 4/' is odd, whereas the amplitude a is even with respect to z. We recall that 

4/V(L^) = -4/V(L^ + 27Tfa) = 

due to spatial symmetry of the component 4/+. Therefore, the continuity boundary conditions 
ao(L^)4/Q(L|u) = a+(L/i)4/'_)_(L|u) = a+(L/U + + 27r/r) 
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yield the boundary values for o: 


(4.63) 


2ao{Lfj,) = a+{Lfj,) = —a+^Lfi + 27r|u). 


On the other hand, 'h" is expressed by the stationary NLS equation (|4.2p . so that 4'" is continuous 
at the vertex points. Therefore, the derivative boundary condition 

C/^(L/x) = C/;(L/r) - C/;(L/i + 27r/i) 

yields the boundary values for the derivative of a: 


(4.64) 


^^(L/i) [2ao(L^) - a\{Ln) - a+(L/i + 27r;u)] = -ao(T/x)4'" (L^u). 


Thanks to the symmetry condition, we are looking for odd U and even a, so that the conditions on 
a at the other vertex point repeat boundary conditions (j4.63p and (j4.64p . 

After the boundary conditions (|4.63p and (j4.64l) are identified, we substitute the product form 
U{z) = a{z)'^'{z) into the eigenvalue problem £+[/ = \U. After multiplying the resulting equation 
by T'(z), we obtain 


A 

dz 


da 


{¥) - =A{^.') «. 


After multiplying this equation by a and integrating by parts, we obtain 


(4.65) 


aa' (^')^ 


where the squared brackets indicate the total jump at the vertex points: 

[/] := - 27r/i) + /o(T/i) - fo{-Lp) + f+{Ln + lirp) - f+{Lfi). 

To compute the total jump explicitly, we use the symmetry on a and T', as well as the boundary 
conditions (j4.63p and (I4.64h . As a result of straightforward computations, we obtain 

aa' (T')^j = 3a+(L//)2T'+(L/i)T" (L^). 


By Lemma HTTl for sufficiently large fj,, we have T" (L^) > 0 and < 0, therefore, aa' (T')^ 

0. The quadratic form p4.65p implies that the corresponding eigenvalue A is positive. This completes 
the proof of the lemma. □ 

Remark 4.11. Lemmas 14.91 and 14.101 prove the assertions of Theorem 11.21 


< 


Remark 4.12. The method of the proof of Lemma [4.101 is inconclusive for the asymmetric standing 
wave given by Lemma 14.31 Indeed, the total jump condition without spatial symmetry of the 
asymmetric standing wave T are given by 


aa 




= ^a+(L/.)2TV(L/r)4/" (L/r) - ^a_(-L/.)24/'_(-L/r)4/" (-L/r), 


where 'k(j_(±L/i) > 0 and 'I'((.(±L/r) > 0 if ^ is sufficiently large. Therefore, the sign of aa' (T')^ 
depends on the balance between a+{Lfi) and a-{—Lfi) relative to 4'+(L//) and 4'_(—L/r). 


5. Numerical approximations of the ground state 


Here we illustrate numerically the construction of the standing waves of the stationary NLS 
equation (11.51) and the corresponding ground state of the constrained minimization problem (11.71) . 
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5.1. Numerical Methods. To compute solutions of the stationary NTS equation ()1.5p for a given 
A, we will use both a Newton’s method (largely the Matlab based program nsoli), as well as the 
Petviashvilli method 1211129]. Rigorous convergence estimates for the Petviashvili method have 
been established in [ 23 | and recently refined in Eoj. Often, we will use the Petviashvilli method to 
initially land on a branch, then continue using the more delicate Newton’s method machinery. The 
basic approach to the Petviashvilli method starts with an initial guess uq, which will generally take 
to be a Gaussian function centered either at the central link or at one of the two loops. Then, for 
a given A < 0, we construct solutions to the stationary NTS equation (11.51) by defining 

Un+l = M[u„]'^(|A| - A)~^{2\Un\‘^Un) 


with 

M[U] = 

2|F|li4 

and 7 = §• We iterate until \un+i — Un\ < then declare such the final state to be a good 

numerical approximation to a fixed point. Once we have constructed a solution that is centered 
either at the central link or at one of the two loops, the Newton solver can be used to continue that 
branch with great accuracy. 

In order to set up our discretization, we approximate the Laplacian operator A using a second- 
order symmetric finite difference stencil with N uniformly space grid points of size h = 2t: jN on 
each of the two loops of length 27r and M grid points on the interior section of length 2L. In order to 
allow an approximation of the 0{h?) order, we first ensure that the discretized Laplacian operator 
is symmetric by taking L = rmi for mGNorl/mGN and choose N such that M = AL/tt G N. 
To enforce the Kirchhoff boundary conditions (II.ip . we take the higher order difference calculations 


u_{-L) 
u'_{—L — 27r) 


—'Uo (—L 2/i) 4'Uo (—L /z) — 3uo(—T) 

’ 

U-{—L — 2h) — 4u-{—L — h) + 3u_(—L) 

2h ’ 

—U-{—L — 27r + 2h) + 4u_(—L — 2Tr + h) — 3u-{—L) 


which allows us to replace uq{—L) = u-{—L) = u-{—L—2tt) everywhere it appears in the symmetric 
difference for the Laplacian. There is a symmetric argument for the other Kirchhoff boundary 
condition in (jl.2l) . 


5.2. Numerical Findings. Using the graph Laplacian approximated to the second order and 
various methods for constructing solutions of the stationary NLS equation (II.5p for a given A < 0, 
we attempt to verify various properties the ground state branches discussed in Theorems ll.ll and ll.21 
We present the results of our various numerical studies in Figures [6l 171 [8l [TOl andfTTl The code 
for computing these are made publicly available at www.unc.edu/~marzuola/mp_graph_code/, 

In Figure [G] we plot the form of the ground state computed using Petviashvili’s method with 
symmetric initial guess localized in the center link for L = 7r/2 and a variety of A values. For values 
of A = —0.01,—0.1,—1.5,—10.0, we observe the computed ground state go from the constant 
solution dni) as in Lemma 13.11 to the positive asymmetric wave as in Lemma 13.41 then in an 
intermediate region to the positive asymmetric wave as in Lemma [4.3l finally, settling on the positive 
symmetric wave as in Lemma 14.71 This confirms the predictions of ground states in Theorems ll.il 
and 11.21 

Figure [7] plots the form of the stationary solution computed using Petviashvili’s method with a 
loop centered initial iterate for L = 7r/2 and A values —0.01,—0.1,—1.5,—10.0. Compared to the 
outcome on Figure [6l we observe that the positive asymmetric wave remains to exist at least up to 
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Figure 6. Solutions of the stationary NLS equation (ll.Sp for L = 7r/2, when the 
initial iterate is localized on the central link, for A = —0.01 (top left), A = —0.1 
(top right), A = —1.5 (bottom left), and A = —10.0 (bottom right). The values at 
±L are marked with a red circle. 




Figure 7. The same as in Figure [6] but for the initial iterate being localized on the 
right ring. 
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Figure 8. Comparison of the standing waves (solid line) localized in the central 
link (left) and in one of the rings (right) to the rescaled solitary wave profile (dots) 
for L = 7r/2 (top) and L = 27r (bottom) for A = —10.0. 




Figure 9. The bifurcation diagram for two positive solutions localized in the central 
link and in one of the two rings as well as for the constant solution for L = 2tt (left) 
and L = 'K!2 (right). 

A = —10. Thus, both positive waves (|1.9I) and (ll.lOp coexist for large negative A as in Theorem 

o 

Figure [8] plots standing waves for A = —10 localized in the central link and in one of the two rings 
and compares them to the appropriately rescaled solitary wave (|4.3I) . The agreement illustrate the 
representations (|1.9p and (jl.lOl) in Theorem 11.21 

Figure [9] shows the evolution of the charge Qo with respect to A for the constant solution and 
the two positive waves continuing from the constant state for L = 2tt (left) and L = 7r/2 (right). 
We compute these branches by using both Petviashvili’s method for sufficiently large |A|, then 
continuing with Newton’s method towards the constant branch. The bifurcation of positive waves 
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Figure 10. Left; comparison of the bifurcation diagram for positive solutions lo¬ 
calized in the central link and in one of the two rings for L = 2tt. Right: the second 
eigenvalue of the linearization operator L+ associated with the two positive solutions 
plotted with respect to A in a semilogy plot. 




Figure 11. The same as Figure [TO] but for L = 7r/2. 


from the constant state is the pitchfork bifurcation as shown in Lemma [3.41 and claimed in Remark 

ESI 

In Figures [Till and fTTl we numerically continue to large values of |A| and explore two properties 
discussed in Lemmas 14.91 and 14.101 for L = 2tt and L = 7r/2 respectively. The left panels show 
the bifurcation diagram zoomed for an interval of large |A|. The two positive waves are computed 
via Newton iteration from peaked states at A = 10.0 concentrated in either the loop or the center 
link. In agreement with Lemma 14.91 we can see that the positive symmetric wave becomes the 
ground state of the constrained minimization problem (USD for a large negative A, because it has 
the smaller value of Qq for a fixed A. We also observe that the values of Qq are much closer for 
large L than for small L. Although only a short interval of A values is shown, the same trend 
continues all way to A = —10. 

The right panels of Fieures [TOl and 1 11 1 show the evolution of the second eigenvalue of the linearized 
operator L_|_ for both the positive waves. In agreement with Lemma 14.101 the positivity of the 
second eigenvalue of L+ for the positive symmetric wave represented by (|1.10l) is observed to be 
quite robust for both L = 2tt and L = 7r/2. However, we also observe that the second eigenvalue of 
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L_|_ for the positive asymmetric wave represented by (|1.9p is also positive, which implies that the 
positive asymmetric wave is also a local constrained minimizer of energy at a fixed charge. 


Appendix A. Proof of Proposition 14.61 


Proposition 14.61 generalizes formal asymptotic expansions given by formulas 16.15 in [T]: 

(A.l) sn(^; A:) = tanh(^) + ^(1 — k^) [sinh(^) cosh(^) — sech^(.^) + 0((1 — 

(A.2) cn(^; A:) = sech(,f) — ^(1 — A:^) [sinh(^) cosh(,^) — tanh(,^)sech(^) + C>((1 — A;^)^), 

(A.3) dn(^; A:) = sech(^) + -(1 — k^) [sinh(^) cosh(^) + tanh(^)sech(^) + 0((1 — A:^)^), 

where the expansion is understood in the sense of the power series in (1 — k‘^) as A —>• 1 uniformly in 
In the same sense, formulas (|4.40p . ()4.4ip . and (j4.42p follow from expansions (jA.ip . (|A.2p . and 
()A.3p . Here we give a rigorous proof of these asymptotical representations, as well as the bound on 
the first derivative given by (I4.43|) . We only give the proof for the Jacobi elliptic function dn. The 
proof for other Jacobi elliptic functions is similar. 

First, let us recall the basic identities for the Jacobi elliptic functions 

(A.4) dn^(^; A) + A^sn^(^; A) = 1, cn^(^; A) + sn^(^; A) = 1, 


from which it follows that dn(^; A) is monotonically decreasing from dn(0; A) = 1 to dn(iF(A); A) = 
Vl — A^, when ^ changes from 0 to K{k), where K{k) is the complete elliptic integral. Also recall 
that dn(^; A) is an even, 2iF(A)-periodic function of ^ for every A G (0,1). 

From the integral representation 8.144 in m, for every A G (0,1) and every ^ G (0,iF(A)), we 
have 


(A.5) 


e = 


dt 


/dn(5;fc) Vl - - 1 + A^ 

Let uo(0 := dn(^;l). Using the parametrization t = sech(x), we obtain from (IA.5P for every 
^ G M+: 

dt r 

Jo 


fXO 




dx = xq, 


JuoiO iVl - V 

where xq is a positive root of sech(x) = uq. Therefore, uo{^) = sech(xo) = sech(^) as in the first 
formula of (j4.42p . 

Differentiating (IA.5I1 in A and using the identities (|A.4p . we obtain 


(A.6) 


9fcdn(^; A) = —A^sn 


(VA)cn(^;A) [ 

J dn 


dt 


dn(^;fc) Vl - - 1 + A^)^ 


Let uo(0 := ^fcdnP^;!). Using the first formulas in (I4.40p and (I4.4ip . as well as the same 
parametrization t = sech(x), we obtain 


(A.7) 




tanh(,f)sech(^) / 

J lin 

tanh(,f)sech(^) 

Jo 


UQ 


dt 

(i) VVi - V 
cosh^ {x)dx 


— i tanh(^)sech(^) [sinh(^) cosh(^) + , 


which justifies the second formula in (|4.42l) . 
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It remains to justify the bound (14.431) on the first derivative of u{^-,k) := dn(^; fc) in k for 
^ E (OjVr^u) and k E {k^{iJ.),l), where /r is sufficiently large. We recall that u{-;k) satisfies the 
second-order differential equation for every A; E (0,1): 


(A.8) 


cP 


u 


- ^ + (2 - k'^)u - 2u^ = 0 , 


CG{-K{k),K{k)). 


Let us introduce the linearization operator 

(A.9) Lk := + 2-k^- Qu\i-k) : H^,,{-K{k), K{k)) ^ K{k), K{k)). 

As k ^ 1, the linearization operator defined by (IA.9P converges in some sense to the limiting 
linearization operator 

(P 

(A.IO) Lk=i := + 1 - 6sech2(e) : ^^(M) ^ L\R), 

which has a negative eigenvalue at —3 associated with the even eigenfunction Uq, the zero eigenvalue 
associated with the odd eigenfunctions Uq, and the essential spectrum at [l,oo), where no(C) = 
dn(^; 1) = sech(,f). 

For convenience, we drop the first argument ^ in the definition of u{p, k) = u{k), so that we can 
denote the partial derivative of u with respect to A; by a prime. Then, we note that 


(A.11) Lku\k) = 2ku{k). 

Since d^u = —fe^sn(^; A:)cn(^; k) solves L^d^u = 0, we obtain the unique 2Ar(fe)-periodic and even 
solution of the differential equation (jA.lip by variation of a constant: 

(A.12) u'{k) =-ksn{p,k)cn{p,k) [ 

Jo cnP^'-,k) 

This representation is complementary to (IA.6P and it admits the same expression as in ()A.7h if the 
limiting values of the Jacobi elliptic functions in (lA.ll) . (|A.2p . and (|A.3p are used. 

Now we note that although the inhomogeneous equation (|A.11I) can be uniquely solved in the 
limit A: —?■ 1 in H‘^(R) because uq is orthogonal to ker(Lfc=i) = span(uQ), the limiting expression 
for u'{l) contains an exponentially growing function of ^ as per the explicit expression ()A.7p . This 
is because the homogeneous equation Lk=iv = 0 admits an even exponentially growing solution 
v{p) = cosh(^) J- v{p), where v E FA^(M). 

Using further differentiation of (lA.llI) in A:, we obtain a chain of linear inhomogeneous equations 
Lku"{k) = 4A:u'(A:) -|- 2u{k) + 12u{k){u'{k))'^, 

Lku"' {k) = Qku"{k) + 6u'{k) + 36u{k)u' {k)u'' (k) + 12{u' (k))^, 

Lku''"{k) = 8ku"'{k) + 12u"(A:) -|- A8u{k)u'{k)u'''{k) + 72{u'{k))‘^u''(k) + 36u{k){u''(k))'^. 


Inspecting the right-hand sides of these linear inhomogeneous equations in the limit A: —)• 1 and 
inverting Lfc=i on the even right-hand sides, we obtain that all derivatives of u in A: are exponential 
growing functions of ^ with the following growth rates: 

u''(l) = 0(^sinh(|)), u'"(l) = 0(cosh(3^)), sinh(3^)), as | ^ ±oo. 

It follows by induction (the proof is omitted) that for every A: E N, we have 

(A.13) tt^^^’''^^(l) = C)(cosh((2A; J- 1)^)), = C)(^ sinh((2A; J- 1)^)), as ^ Too, 
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where the implicit constants grow polynomially in k. The A^-th partial sum of the Taylor series 

^ 1 
n\ 

n=l 

converges for every ^ G (—7r/i,7r^) and k G where /c*(/r) = 1 — as 

fi ^ oo. Therefore, u'{k) — n'(l) = limAr^oo Sn is well-defined by the majorant power series in the 
corresponding domain. Moreover, if /x is sufficiently large, then 

|n'(/c) - < CjA: - 1|(1-h l^l)cosh(^), ^ G (-tt/t,tt/t), k e {k^{n),l), 

where C is a positive /r-independent constant C. This bound is equivalent to the third bound in 

(IM) . 
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